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The propagation of a narrow-band signal radiated by a point source in a randomly layered absorb- 
ing medium is studied asymptotically in the weak-scattering limit. It is shown that in a disordered 
stratified medium that is homogeneous on average a pulse is channelled along the layers in a narrow 
strip in the vicinity of the source. The space-time distribution of the pulse energy is calculated. Far 
from the source, the shape of wave packets is universal and independent of the frequency spectrum 
of the radiated signal. Strong localization effects manifest themselves also as a low-decaying tail 
^.^ . of the pulse and a strong time delay in the direction of stratification. The frequency-momentum 

(~| ' correlation function in a one-dimensional random medium is calculated. 

a : 

I PACS numbers: 05.40.-a, 05.60.-k, 41.20.J, 84.40.-x 

C/3 ' 



I. INTRODUCTION 



wj-i ' The propagation of quantum wave packets and pulses of electromagnetic radiation in disordered media is a classical 

r^ : problem with a long-standing history. The continued interest of physicists in this problem is stimulated both by 

Q the quest to better understand such fundamental problems of disorder as correlations in momentum-energy space, 

O , localization of time-dependent fields, Wigner time delay, etc., and also by the growing number of applications that 

pulsed signals find in modern electronics, telecommunications, optics, and geophysics. Considerable theoretical and 

^__l \ experimental investigations have been expended to study the propagation of pulses in randomly inhomogeneous media 

J> • in diffusive regime (see, for example, Ref. |l| and references therein). Much less studied is the space-time evolution of 

^^ ' wave packets in disordered one-dimensional and layered systems where the interference of multiple scattered fields is 

OO of crucial importance. 

It was shown in Refs. ^^ that in a homogeneous on average, randomly layered medium where the refractive index 

(potential) is a random function of one coordinate only waves (quantum particles) are localized in the direction 

of stratification and propagate along layers forming the so called fluctuational waveguide. The statistics of wave 

fields radiated by a monochromatic point-like source in a randomly layered medium was studied in Refs. 0,0. For 

"^ I its analysis, the resonance expansion method was applied to calculate correlation functions of plane harmonics with 

d ' different "transverse energies" , i.e. squared projections of the wave vector on the axis of stratification. In the case 

S , of a non-stationary signal, the problem becomes much more complicated because it involves correlation analysis of 

I ' waves with different both frequencies and transverse wave numbers. 

' ^ , In the present paper, we investigate analytically the space-time distribution of the average intensity of pulse field 

^ ' radiated by a point source in a randomly stratified weakly scattering medium with dissipation. As an intermediate 

fj \ result, the correlation function of the propagators (Green functions) with different frequencies and transversal wave 

• • ' numbers is calculated. Localization of the constituent plane harmonics is shown to result in channelling of the pulse 

. ^ ' within the fluctuation waveguide and in a significant modification of the spectral content of the signal far away from 

^ . the source. The shape (envelop) of the pulse in the far zone is calculated. It is shown to be universal and independent 

V^ ' on the spectrum of the radiated packet. This is due to both the filtration of the harmonics by their localization 

radii transverse to the layers and to the difference in phase velocities of those harmonics in this direction. The same 

reasons cause strong time delay of the pulse when the receiver is shifted towards the direction of stratification from 

the horizontal plaiie in which the source is located. This effect is a clear manifestation of the delay time concept 

introduced earlieiu on the basis of scattering phases of quantum particles moving in disordered media. 

II. FORMULATION OF THE PROBLEM 

We consider the wave equation for the scalar non-monochromatic field G(R, Ro|t) radiated by a source located at 
a point Ro in an infinite medium which is randomly stratified along z-axis, 
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G(R, Ro|i) = 47rJ(R - Ro)A(t)e 
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(1) 



Here A is Laplacian, e{z) = Eq + 5£{z) is the (random) dielectric permeability with the mean value Eqi f is the 
conductivity of the medium, A{t) is the envelope of a wave packet (pulse) with the carrier frequency loq. In what 
follows we consider a narrow-band wave packet which means that A{t) is a smooth function as compared to the 
oscillating exponential in the rhs of equation (|l|) . 

Since permeability e in Eq. (nl) depends on one coordinate only, the problem of finding the mean intensity (/) = 
\\G\ ) at a given point R is reduced, after Fourier transformation in {x, j/} plane, to calculation of the one-dimensional 
correlator of harmonics propagating along the axis of stratification. 
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Here, angular brackets denote statistical averaging over the ensemble of random functions Se{z), Jq (x) is the Bessel 
function, g = Rn — Rqii is the in-plane distance from the source, Rii is the radius-vector component parallel to the 
layers, f{uj) is the spectral function of the radiated pulse. Function K{. . . ) in Eq. (g) is the two-point correlation 
function 

K (z, zo|ki, K2; ^1, ^^2) = (G (z, zo\ki,uji) g* (z, zo|k2, ^2)) (3) 

that we present in a form more convenient for subsequent calculation by changing the integration variables, viz. 



K{z,zo\ki,K2;uj,^) = {Q {z,ZQ\q^ + lS.q^,uj + n) Q* {z , za\q^ , uj)) . 
Here uo — llI2i f^ = cji — lj2, q^ = k^^ — k\, k^j = eo (^/c) • The "energy" difference /S.q^ is given by 
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7(0;) — i-na /uj. In Eq. (jj), Qiz, ZQ\q'^ ,uj) is the Fourier transform over in-plane coordinate R|| — Ro|| and time of the 
Green function G from Eq. d|). This function obeys the equation 



-— + q^ + iQ + 5e [z] - 



dz 



t; (z,zo|g^,cj) = 47r(5(z- zo) , (6) 

Formula (|g) for energy difference Ag^ is valid, strictly speaking, in the case of weakly dissipative medium, i.e. when 

|7(^)| « 1 . (7) 

Under the assumptions of weak dissipation and spectral narrowness of the pulse the expression (0) can be readily 
transformed to the form 



Aq^ ^ kI- kI + 2£q- 



LjjVL 



-27H(^)^ 



(8) 



Correlation functions of the type (H), (0) with uji = lu2 (i.e. with fJ = 0) that aiipear in the theory of stationary 
processes in ID disordered systernSj^an be calculated using diagrammatic methods,LlLl functional method of Ref. ||, or 
the resonance expansion method.ou In the next section the latter approach is shown to be quite universal and well 
applicable (after some modification) also to non-stationary stochastic problems, in particular, for calculation of the 
correlation function (0) with ft ^ 0. 



III. RESONANT SCATTERING APPROXIMATION FOR FIELD CORRELATORS 

In this section, the method used in Refs. pP for calculating statistical moments of the field radiated by a monochro- 
matic point-like source is generalized to the case of pulse signals. The method allows for rigorous calculation of the 
correlator (H) provided a single scattering can be regarded as weak. 

A. The resonance expansion method 

To calculate the intensity using Eq. (^) we have to know the Green functions in (||) for all values of q^ in the 
interval — oo < q^ < -I-cxd. However, as it was shown in Ref. m in the case of weak scattering the contribution of 
spatial modes with q^ < (so-called evanescent modes) is significant only in a thickness of |z — zg] ~ k~^ near the 
source position zq. In the rest of space the intensity is largely determined by the propagating (extended) modes for 
which the Green function obeys equation (g) with q^ > 0. The key point of the following calculations is the so-called 
resonance expansion of this Green function, 

e(z,zo|'?',o;) =ei(z,zo)e^^("-"")+e2(^,^o)e-^''("-"")+e3(;2,^o)e^«^^+^'')+e4(^,^o)e-'''("+^°) , (9) 

where Gi(z,zo) = Qi [z , ZQ\q'^ , oj) are smooth factors in comparison with the "fast" exponentials. The assumption 
of smoothness of the "amplitudes" Qi{z,zo) is based on the requirement for weak scattering (WS) of the pulse- 
constituting plane harmonics, which means that the extinction lengths L of the harmonics, see Eqs. (|l4| ) below, are 
large compared to their wavelengths and to the correlation radius Tc of 6e{z) as well. 

Formula ()^ represents the exact Green function as a sum of relatively small packets of spatial harmonics centered 
at four basic ones, viz. exp[±ig(z ± zq)]. Such a form of the solution of Eq. (o ) implies that only resonant harmonics 
in the power spectrum of the permeability fluctuation 5s{z) contribute significantly to the scattering of a wave with 
the wave number g, namely the harmonics with the momenta close to zero, which are responsible for the forward 
scattering, and close to ±2q (backward scattering). 

Using Green function in the form (pi) one can perform spatial averaging of equation (pf) over an interval 21, such 
that q~^, Tc ^^21 <^ L. As the result, for the matrix 



of the smooth amplitudes from IM) the equation follows 
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dz 
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Here a^ and a are 2x2 matrices 
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superscript (f) stands for Hermitian and the asterisk for complex conjugation, respectively. Random functions ("po- 
tentials") ri{z) and C(z) are constructed of narrow packets of spatial harmonics of Se{z) as follows. 

On the assumption of weak scattering, functions 77(2;) and Ciz) can be thought of as Gaussian random processes 
irrespective of the statistics of Se{z)B Correlation of these functions was studied in detail in Ref. |l^ where the evidence 
was given that only two binary correlators of the potentials dl3), viz. (77(2)77(2')) and (C(z)C*(z')), are different from 
zero, and can be replaced by weighted (5- functions, 

(7;(z)77(z')) = Lj'Siz - z') , (C(^)r (^')> = L^'Siz - z') . (13) 



In (|lj), length parameters Lf^t, are given by 



(2g)2 ^ ^ ^ /cx4 (2q)2 
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W{p) is the Fourier transform of the binary correlation function of the permeability fluctuations, 

W (z - z') ^ {6e{z)Seiz')) . (15) 

It is shown in Ref. O that (n^ are nothing but the extinction lengths related to the forward (/) and backward (6) 
scattering of the harmonics with the wave number q and frequency uj. In terms of these lengths, the WS conditions 
used when deriving equation (O) are expressed through the inequalities 

g-\rc<L/,6. (16) 

Notice that the relationship between small lengths q^^ and Vc is of no crucial importance. It only specifies the Fourier 
components of the correlation function ( |l5| ) and thus the possible distinction between the "forward" and "backward" 
extinction lengths dl4). 

For the resonance approximation (equivalent to WS requirement) to be most efficient it is advantageous to represent 
both of the Green functions entering the correlator (|j) in the form of the expansion (BJ) with the same fast exponentials, 
i.e. with the same wave number q. Although the amplitude functions Gi in (H) cannot be found explicitly, this 
representation proves to be quite helpful for the calculation of the correlator (|j). Indeed, if we present both of Green 
functions from (Q) in the form of expansion (^) with the same fast exponentials, only "diagonal elements" of the 
product GG* remain non-zero after the averaging (see next subsection). 

Note that functions Gi in (g) can be recognized as slowly varying amplitudes if along with the WS condition ( p^ ) 
the inequality holds 

|Aq2|<g2. (17) 

Physically, this inequality is natural for the definition of weak scattering since the quantity q^ has the meaning of 
energy in Eq. (||). As it will be clear from the subsequent calculation, the inequality ( [l7|) is coincident with the 
conditions (ttq) supplemented by the requirements of weak dissipation, 7(^0) ^ Ij and narrowness of the pulse 
frequency band. 

Green function of equation (0) is the solution of a two-point boundary-value problem with conditions given at 
z — > ±cx). However, to systematically perform the averaging over random potentials without resort to finite-order 
perturbation approximations (that fail to take into account correctly the interference of multiply scattered waves in 
one-dimensional random systems) it is much more convenient to deal with random functions obeying Cauchy problems 
which are causal functionals of the random potentials (^2|). Fortunately, the elements of Green matrix (|l^) can be 
factorized into products of the auxiliary one-coordinate functions, each meeting the initial- value problem conditioned 
at either plus or minus infinity. The factorization scheme is outlined in Appendix |A|. Evolutional character of the 
equations for those functions allows to obtain finite-difference equations ( pi| ) (see Refs. p^p^ ) for auxiliary correlators 
with the help of which the correlation function (jj) can be appropriately calculated. 

B. Asymptotics of the correlation function Eq. (W) 

To obtain analytic expressions for the correlation function (0), we assume the medium to be statistically uniform 
on average in z direction, and then pass from the coordinate representation of (jj) to its Fourier transform over the 
variable z — zq, 
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After substitution of the matrix elements of Green matrix (hOf) in the form ( A7 ) into (W) , and then into (p_8|) (where all 



non-coordinate arguments of function (Eh are omitted for a while), one has to expand the factors A'^{zo) and A*{zo) of 
those elements in series of the products r^r_ and F^Fl. Possibility of such expansion is ensured by retro-attenuation 
iO in Eq. (p). In the course of statistical averaging, each of the terms of the double series produced for the function 



K{s) is decomposed into a product of functionals of different causal types, viz. "plus" and "minus" type. Since the 
potentials ( jl^ ) are effectively (5-correlated, the supports of the random functions entering the functionals of different 
types do not overlap. Therefore, statistical averaging of those functionals can be performed independently. 

It can be shown that the correlators of th e ty pe ((F^) (Fj.) ) with n ^ to in the double series for K{s) are 
exactly equal to zero. Indeed, from equations (A3) it follows that the functional series for the functions 7r-|-(z) consist 
solely of the terms with equal numbers of the functional factors C and C*: whereas the quantities 7± contain extra 
factors, C* for 7+ and (* for 7_. Inasmuch as under WS conditions dlq) functional variables (03) can be regarded as 
Gaussian-distributed random fields, the above-indicated correlators have non-zero values only if n = to,. 

The foregoing procedure has been described in detail in Refs. 12,13. Omitting here tedious calculations we present 
the final result of the averaging. The function K{s) is represented as a series. 
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where i?„ and P„(s) are the auxiliary correlation functions of the form 



i?„ = ([r^(z)r;(z)] 
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that obey the following finite-difference equations 

(3Rn — n {Rn+l + Rn-i — 2i?„) 
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(20b) 
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Pn{s) - Pn-l{s) + ( isCb + ^ ] Pn{s) + l3nPn{s) 



2n' + 2n + 1 + -^(1 + 2n)^ P„(s) = A (i?„ + i?„+i) , 
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where the following notations arc used 
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Equations ( pi| ) have to be supplemented with the re quire ment that functions i?„ and Pn{s) tend to zero as n — > 00. 
As regards their behaviour at n = 0, from definition (|20a| ) it follows that Rq ~ 1 whereas for Po{s) integrability over 
the variable s is sufficient. Equations of this type have been studied in Refs. ||J3J§ and p^ljl^ in the context of the 
conductivity of ID disordered systems. 

The terms proportional to a^ allow, in principle, for arbitrary non-stationarity of the wave to be taken into account. 
Yet narrowness of the pulse frequency band assumed in this paper is consistent with the inequality |a;| ^ 1 allowing for 
equations ( |2l| ) to be solved perturbatively in this parameter. In Appendix B, it is demonstrated that if the inequality 
holds 1/3 1 ^ 1, the summands oc a^ in Eqs. (|2^) contribute negligibly to the sum (p^, in accordance with the condition 
(p^). As a consequence, in the limit of |/?| — + we arrive at the following expression for correlation function (0), 



K{z,zo;ki,k2;uj,Q) 
Here the notations are used 
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2 cosh^ (7r/x/2) ' 

In the limiting case \P\ ^ 1, the terms proportional to a^ in Eqs. ( |2l| ) result in small, though not a priori negligible 
corrections to the basic approximation for the correlation function H), 
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It will be shown in the next section that the average intensity of a narrow-band signal is mainly determined by the 
behaviour of the correlator (0) at ki ~ K2, that corresponds to |/3| <C 1 and, consequently, to the asymptotic expression 



IV. CALCULATION OF THE PULSE SHAPE 

To calculate the average intensity (/ (R, Ro|i)) we evaluate the integrals over fi and lu in Eq. (ph with the correlator 
K (z, zo|ki, K2', oj, ^) given by formula (E3) and function ip{lli + fl — ojq) presented in the form 



ip{uj + n- ujo) = / dt'A{t') eK-'+f2--'o)t' _ (26) 

From asymptotic expression ( p3|) it follows that in the lower half-plane of the complex fl their is a pole with Im Up — 
—Ania/eo, {Aq = 0) which allows to calculate the integral over f2 as a residue in this point. Is is obvious that a 
non-zero result is obtained for those values of t' that are limited by the condition 



t' < t*{q2,uj) = t-T[q2,uj) 



(27) 



Here c = c/y^, notation T{q,Lo) stands for the time interval necessary for the plane harmonics q to pass from the 
plane of the source (2:0) to the plane of the receiver (z), 
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The next step is to calculate the integral over w in Eq. (g). Due to the presence of the narrow function (^*(w — wq), 
all physical quantities in the integrand, in particular 7(0^) and Km(a;), can be taken at the carrier frequency wq- Since 
in the present paper we are mainly interested in localization effects, dissipation in the medium is supposed to be small 
enough, and the dissipation rate of the carrier harmonics is much larger than the corresponding extinction lengths 

(0), 
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Subject to the condition (|29|), the integral over u) recovers the function A*(t'), whereupon the average intensity is 
reduced to 
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From here on we address the case when the receiver is located far from the source, so that 

KynQ > 1 . 



(31) 



To integrate in Eq. (|30| ) over qi2 we use the integral representation of the Bessel functions and implement the 
saddle-point method. Both of g-integrals in (BOl) have the same saddle point 
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so that simple calculation then yields 
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(33) 



This result for the space-time distribution of the average intensity of a point-source-radiated narrow-band signal is 
rather general and is valid for arbitrary envelope A{t). It is simplified substantially when the distance g is large 
enough for the pulse duration to be less than the time of the pulse arrival at the observation point in homogeneous 
{5£{z) = 0) medium. In this case the upper limit in the integral over t' in (|30| ) can be extended to the infinity, all 
functions in the integrand of (p3[) can be taken at t' = 0, and from (p3) we obtain 
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Here Lm = Lb{Km,i^o) is the largest value of the backscattering-induced extinction length (localization length) 
corresponding to the most "energetical" (i.e. q2 = «,„) harmonics, and 



T: 



dt' \A{t')\^' 



Although the intensity of a monochromatic field is known to be a strongly fluctuating, not a self-averaged quantity in 
ID disordered systems, the integration of the correlator (23) over parameters u and q (the last integration corresponds 
physically to the summation of plane harmonics with different angles of propagation) serves as an additional averaging 
factor that suppresses fluctuations of the intensity of the wave packet, and therefore makes the results obtained by 
ansemble averaging, (63) and (pih, more physically meaningful. 



DISCUSSIONS OF THE RESULTS 



Equations (p3[) and (B4J) present the space-time dependence of the average intensity of a narrow-band pulse signal 
radiated by a point source in a randomly layered weakly dissipative medium. Here we dwell on the main physical 
characteristics of the result that are manifestations of the strong Anderson localization in ID disordered media. 

First, it is evident from ( p3| ) and (Q) that the pulse field is exponentially localized in z-direction within a 4Lm 
thick layer whose central plane z = zq is the plane where the source is located. In other words, the point-source pulse 
radiation is channelled, much as the monochromatic radiation is, within the fiuctuation waveguide which is created 
owing to the interference of multiply back-scattered plane harmonics even though the regular refraction does not exist 
in the system. 

The next interesting feature is that the narrow-band pulse acquires universal shape at large distances in the 
fiuctuation waveguide. Indeed, under the assumption that in-plane distance to the observation point is such that 
p > cT, the average intensity is described by formula ( |34| ) and depends on the envelop of the incident pulse, A{t), 
only through the normalization constant T. In Fig. m, a set of curves is presented depicting the intensity distribution 
as a function of the in-plane distance g at a given time t for different distances in the direction perpendicular to 
layers, C = |z — zo\ /4Lm. In Fig. the time dependence of the pulse intensity is shown at a certain distance g 
and three different (. As it is seen from the graphs, during the propagation in the fiuctuation waveguide the signal 
acquires a rather slowly decaying tail, and at large distances from the maximum of the pulse the intensity decreases 
in time proportionally to i^^. The weak sensitivity of the wave packet to its initial shape is due to the fact that in 
randomly layered media (in distinction to free space) a point source radiates-anly those eigen modes that are localized 
in a narrow (of the size of the localization length) stripe near the sourccaa The (random) set of these modes is a 
fingerprint of each realization of random potential and is independent on the way of excitation. 

Another peculiarity of the pulse propagation in a randomly layered medium is a sort of "anisotropy" of the time 
delay of the wave packet: the arrival time increases with increasing \z — zo\ faster than it does when g grows. Indeed, 
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FIG. 1: Instantaneous spatial distributions of the intensity in the fluctuation waveguide at difi'erent ^ ; 
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FIG. 2: Time dependences of the intensity at a given in-plane distance g and different (^ = '^j^^°' 



the earliest signal arrival time at a point {g,z} is of order 
distribution of pulse in the fluctuation waveguide (|z — xqI < L, 
factor 



\/g^+(z-2o) 



At this moment, if g ^ L,„, the spatial 
given by Eq. (pj) contains the exponentially small 
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The moment, <,„, when the signal 
at the point {e,z} reaches its not 
exponentially small maximum can 
be roughly estimated by equating 
the localization exponent in (p3) to 
unity. Ths procedure yields 
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(36) 



with some numerical coefficient S ~ 
1. Although the estimate ( |36| ) can- 
not claim for satisfactory accuracy, 
it gives a reasonable idea of the 
pulse delay in a stratified medium. 
To accurately calculate the depen- 
dence of the arraival time of the 
pulse maximum as a fuction of the 
transverse displacement | z — zq I one 
should use Eq. ([3^). The results of 
this (numerical) calculation is shown 
in Fig. |. 




FIG. 3: Arrival time of the pulse maximum, tm (in units g/c), vs C, ■ 
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VI. 
CONCLUDING REMARKS 

To summarize, in this paper the problem of the space-time distribution of the average intensity of a narrow-band 
pulse which is radiated by a point source in 3D randomly layered medium has been solved by means of the generalized 
resonance expansion method. The pulse field is shown to be localized in the direction of stratification and channelled 
parallel to the layers within the fluctuation waveguide whose symmetry plane goes through the source location. The 
waveguide originates exclusively from the interference of random fields multiply scattered by weak fluctuations, with 
no regular refraction present in the medium. The typical width of the waveguide is of order of the localization length 
of the harmonics with largest allowable momentum along the axis of stratiflcation. Random lamination of the medium 
leads to a substantial distortion of the pulse shape. Specifically, far away from the source the narrow-band pulse of any 
original spectral content, being locked within the fluctuation waveguide, spreads into a signal with the envelope given 
by the universal function (^) depicted in Figs. [^ and ||. In contrast to homogeneous media, the dependance of the 
arrival time of the pulse maximum on coordinates is strongly anisotropic: it increases drastically as the observation 
point moves in z-direction away from source. This delay is not due to the increase of the path length of the signal, as 
it is, for example, in media with regular refraction, but is caused by the multiple random scattering of the saddle-point 



harmonics (32) 



APPENDIX A: CALCULATION OF THE GREEN MATRIX EQ. (0) 

To find the ID Green function of Eq. ra) we express it via solutions of the appropriate Cauchy problems, 

giz,zo) ^W-' [iP+{z)i;^{zo)9{z - zo) +ij+izo)ib-{z)e{zo - z)] . (Al) 

Here, the functions tp±{z) are the linear-independent solutions of homogeneous equation (o) with boundary conditions 
given at either "plus" or "minus" infinity, depending on the "sign" index, W is the Wronskian of those functions, 6{z) 
is the Heaviside unit-step function. 

In the case of real q, it is reasonable to represent the functions ^±{z) as superpositions of modulated harmonic 
waves propagating in opposite directions of z-axis, 

4'±{z) = 7r_|_ (z) exp{±iqz) — 17^ (z) exp{'^iqz) . (A-2) 

The upper index A indicates that the corresponding functions is rel ated to the first Green function in Eq. 



Under WS conditions, (16), the "amplitudes" 77^(2) and J±{z) in ( |A^ ) are smooth functions in comparison with the 
nearby standing fast exponentials. By averaging the equations for tp±{z) over the interval of the length 21 intermediate 
between "small" and "large" lengths of Eq. (O) we arrive at a following set of dynamic equations, 

±^f (z) +z[v{z) - Aqy2q]7:^{z) + C±(^)7± (^) = , (A3a) 

±7^'(z) -t[r^{z) - Aqy2q]j^{z) + a(^)7r^(^) = . (A3b) 



The function r]{z) in (A3) coincides with the analogous function from (p3), with lu being replaced hy to + fl. The 



functions C±(^) s-^e given by 



Sommerfeld's radiation conditions at z -^ ±00 are reformulated as the "initial" conditions for the smooth amplitudes, 

lim t:^{z) = 1 , lim 7^(2) = . (A5) 

z — >±oo z — )-±oo 

In a similar way the second Green function in ^ is to be represented, keeping in mind that in this case Aq^ ~ and 
= 0. 



Wronskian W in (Al) within the WS limit reduces to 

W = 2iq[n^{z)n^{z)+j^{z)j^{z)] . (A6) 
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By inserting then (A2) and (A6) into ( |Al[ ) and comparing the result with Eq. (||) we obtain for the matrix elements 
ofdlT^ 






g3A{^,Zo\uJ + n,q^ + Aq^)= -—A^{zo) 



^±^ 



0± A 



9^rt(zo)^ 









±^^^r^(zo) + e^^ 



Here 9± — 9\^± [z — zq)] and the rest of notations are 



(A7a) 
(A7b) 






(A8a) 
(A8b) 



The upper sign indices in (A7) correspond to Qi and f/3 whereas the lower signs to Q2 an d Q j, respectively. The 
functions r^(z) and 7r^(z) obey the Riccati-type coupled equations resulting directly from (A3), 



± — i^ = 2i 



± 



dz 
d 1 

rfi7r^(z) 



77(z) 






21 



r^(z)-CiW + C±W[r^W]' 






1 ^^r^^±(") 



(A9a) 
(A9b) 



The initial conditions for Eq s. ( | Ag| ) follow from ([A5|). 

Transition from equations ( |A3D to (A9) is motivated by the following. In the stationary and non-dissipativc limiting 
case (Ag^ = and 7 = 0) the functions r±(z) represent the amplitude reflection coefficients for the harmonics q 
incident on the ID disordered half-spaces (z,±oo), respectively. In the presence of arb itrary weak dissipation these 
functions become modulo less than unity allowing for the factor A^{z) given by Eq. (A8a) to be expanded into a 
series in powers of the product F^F^. Averaging then termwise the double series into which the product of the Green 
functions in (|j) is expanded we arrive eventually at the expression (|l9|). 



APPENDIX B: ANALYSIS OF THE FORWARD SCATTERING CONTRIBUTION 



Equation (21a) can be solved rigorously at a = (see, e.g., Ref. g), therefore it is not difficult to obtain an 
asymptotic expression for i?„ at \a\ <C 1. With the accuracy of the ffist order in a^ one can find that 



R. = p[dte-^^{J^^ 



l + aMl + 2^|< 






l+i-y(3 + 2i) 



(BIO) 



When |/3| <^ 1, the domain corresponding to t ~ |/3| ^ is significant in the integral (1B1C| ). Therefore the contribution 
of the terms proportional to a^ is of the order |a//3p. It follows from Eqs. (||) and (p2D that 



< 



LuCb 



(Bll) 



In rhs of (Bll) there is nothing but a small WS paramet er wh ich governs all the approximations made in the course of 
solution. Therefore, the terms proportional to a^ in Eq. ( [B10| ) lead to the corrections that are less than the calculation 
accuracy, and must be omit ted in (BIO). 

To analyze the equation ( |21b| ) we present the correlation function (|lj), using ( BIC ) at a = 0, in the integral form. 



27r\^ p/3 /•"o 



(B12) 



Here 2/^(0 i^ the function related to the generating function Ys{(,) = J2'^=o ^^Pn{s) by the equality 
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..(0=ln(l-l 



From Eq. ( 21b ) it follows that i/s(0 obeys the differential equation 



(B13) 



di^ d^ ^ ^ d^^ ^ ^ d^\d^ 



l]+i[s 



~2^ 



ysiO + l(^f] QpVsiO 



= A 



2 + (2^ - /3)e«Ei(-0 + ( ^ ) PpiO 



(B14) 



where 



PpiO 



l + 2|^]e^ 



^ i' + ii-P) 



K'-(C'-/3)(2e' + /5)] 



(B15) 



and the differential operator Qp has the form 



Q^- ^l^^) -^ 4^'{|'^J, 






+ ^^^d?^'^ 



^(4^) -^H4^A4'4' 




^^i^-v 



(B16) 



The brackets {•••,•••}+ in ( B16 ) denote an anticommutator. It is evident from ( B12 ) that the solution of Eq. ( B14 ) 
is of importance in the domain C ^ 1 • In that region the estimation is valid 



\\Qp\\-\ppm-^ 






Thus we conclude that the terms proportional to Q^ and P/^(0 , that contain the forward-sc atter ing parameter Cf, 
contribute negligibly, in accordance with WS parameter ( Bll ), to the solution of equation (B14), just in the same 
way as the terms oc a^ in (BIO). 
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